Unfactorizing Polychromatic Penguins. 



Santi Pcris a 

a Grup de Fisica Teorica and IFAE, Universitat Autonoma de Barcelona, 08193 Barcelona, Spain 

Firstly, I report on a recent work -carried out with T. Hambye and E. de Rafael- on how to use the large- 
N c expansion to calculate unfactorized contributions from the strong penguin operators, and their impact on 
observables such as e'/e and the AJ =1/2 rule. Secondly, based on work done with M. Golterman, I explain 
how this calculation allows one to predict some rather dramatic consequences for quenched QCD. This may help 
explain the present discrepancy between lattice and experimental results for e'/e. The emphasis of this article is 
put on the explanation of the method of calculation used, which is fully analytic. This allows one to build some 
intuition and understand the role played by the different hadronic scales in determining the size of the different 
contributions. 



Due to the disparity of scales between the kaon 
and the W masses, the decay process K — > itit has 
very large logarithms, such as log Mw/Mk- The 
technique of Effective Field Theory allows one to 
resum all these large logarithms as the renormal- 
ization group running of certain new couplings ( 
the so-called Wilson coefficients). These coeffi- 
cients appear in the Effective Lagrangian where 
the W and all the other heavy particles of the 
Standard Model are no longer present, i.e. they 
have been "integrated out". Therefore, in this 
Effective Lagrangian only the u,d, s quarks and 
the gluon fields appear explicitly, the effect of the 
heavy particles being encoded in a set of 10 Wil- 
son coefficients and equal number of four-quark 
operators 1 , i.e. 

Ceff = J2 C ^) ^ B {H) , (1) 
AB 

where Qab(m) = qT A q(x) qT B q(x) and Ta,b 
are matrices in Dirac space (color indices are 
suppressed). This Effective Lagrangian, C e ff 
changes the short-distance properties of Green's 
functions with respect to the full Standard Model, 
where the heavy particles were explicit fields in 
the Lagrangian. A set of matching conditions en- 
sure that the physics does not change, however. 
These matching conditions define Wilson coeffi- 
cients and operators and make the whole con- 

1 See, e.g., ref. for the full list of AS = 1 operators. 



struction meaningful so that the scheme depen- 
dence appearing in the Wilson coefficients can- 
cels that appearing in the four-quark operators' 
matrix elements. 

In perturbation theory all this is very well un- 
derstood. Two groups have been able to compute 
the Effective Lagrangian even to two loops 2 . 
However, perturbation theory may be a sensi- 
ble approximation down to the integration of the 
charm quark, but certainly not to describe the 
matrix elements of the kaon. In fact, since the 
kaon is a (pseudo) Goldstone boson, its dynamics 
is appropriately described in terms of yet another 
Effective Field Theory. In the AS = 1 sector, this 
is given by the Chiral Lagrangian 2 : 

gs£& + 527/327 + e 2 g ew Tv (uxf^rfQ^] , 
where 

i=l,2,3 

£27 = l&vhi (C) 13 + (£ M ) 23 0C")n , (4) 
with 

= -iF% Ujx^D „U(x) , A^ 32) = 5 i3 S j2 (5) 

2 In the large- N c expansion one has to include the 77' in a 
nonet. 
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and 



Ql = Qr = Q = diag.(2/3, -1/3, -1/3) 



(6) 



The pion decay coupling constant Fq is in the 
chiral limit, i.e. Fq ~ 87 MeV. The matrix 
field U contains the Goldstone fields of the corre- 
sponding spontaneously broken chiral symmetry 
of QCD, and D^U stands for the covariant deriva- 
tive D^U = dfJJ — ir^U + iUl^ in the presence 
of external chiral sources and r M of left- and 
right-handed currents. The Lagrangian J2J can 
be thought of as the one obtained from (JTJ after 
integrating out all the hadrons but the Goldstone 
octet. How does one go about matching the two 
Effective Field Theories Q and ©? 3 At this 
point we find it very convenient to use the large- 
N c expansion to organize the calculation because, 
being a systematic expansion in QCD, it can be 
carried out either with quarks and gluons or with 
mesons as degrees of freedom. So, in a way, the 
large- N c expansion plays the role of a dictionary 
between these two languages. 

In the following I will concentrate on how to 
impose the matching condition for the case of the 
strong penguin operators 

26 = -8 (sLqii) (Sr^l) , and 

q— u,d,s 

Qi = 4 faW) (q L J»d L ) , (7) 

q— u,d,s 

where qL.R — ^(1 T 7s)<7 and sum over color in- 
dices within brackets is understood. Notice that, 
since four-quark operators mix among them- 
selves, it makes no sense to consider the contribu- 
tion from one of them in isolation. However, as it 
turns out, if we are willing to restrict ourselves to 
the leading logarithmic approximation and only 
to those subleading contributions in 1/N C which 
are enhanced by the number of flavors uf, there is 
an important simplification since the Qifi system 
closes and no other operator needs be considered 
through mixing. I will adopt this simplification 
in the following. 

Firstly, let us define the constant gg in the La- 
grangian [21 as the coupling governing the non- 

s Our approach to this problem has also been explained in 

El. 



diagonal "mass term" r2 rf g g a p for any flavor 
q = u,d, s. For dcfmitcncss, we shall take q = u. 
The matching condition stating that the same 
term be obtained with the Lagrangian reads 4 



9%\Q4.,Qe = CeC//) 



F 6 



(8) 



MS 



+ T^%4 I dQ 2 Q 2 W DGRR (Q 2 

{A f°° 1 

i-TTrfei / d Q 2 Q 2 Wllrr(Q 2 ) 
16n F J ° J MS 

The subscript MS reminds one that these inte- 
grals are UV divergent and have to be regular- 
ized and renormalized, using the same scheme as 
for the Wilson coefficients C^e ■ The parameters 
(ipip) and L5 are also renormalized accordingly. 
In Eq. © all the unfactorized contributions, of 
0{np /N c ), are contained in the terms propor- 
tional to the functions Wdgrr and Wllrr- The 
terms proportional to L5 and unity correspond 
to the factorized contribution from Qq and Q4 
-respectively- and, formally, are of 0(N®). 

The functions Wdgrr and Wllrr are defined 
through the connected four-point Green's func- 
tions 



W 



a/3 
DGRR 



(«) = 



d x d y d z e 



iq.x 



(0\T{D, g (x)Gq d (0)R^(y)Bi s (z)}\0) coan . , 



W 



fi fi a. (3 
LLRR 



(9) 



d x d y d z e 1 ' 



(0\T{L^ g (x)Lf(0)R^y)Rl(z)}\0) cwn , , (9) 

after integration over the solid angle in q- 
momcntum space in the manner specified in 
In Eq. © D Sq = s L q R , G qd = q R d L , L^ q = 
SL'Y^qL, R du — dRj a UR, and similarly all the 
others. One then recognizes the pair of fermion 
bilinears which make up the operators Qe^ in Eq. 
@ except that they are located at different space- 
time points. It is the integral over Q in Eq. (JSJ) 
which brings them back to the same point so that, 
in fact, the matching condition © is nothing but 
the statement that in the four-quark Effective La- 
grangian of Eq. Q one can only generate this 
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two-point correlator between r- and r_ by in- 
serting the combination c 6 Q 6 + C4Q4. Since the 
same two-point correlator is proportional to g$ in 
the language of the Chiral Effective Lagrangian, 
this explains the origin of the condition iJSJ. 

It is clear that if one knew the functions 
Wdg(ll)rb. one immediately would be able to 
compute gs through Eq. (JSJ. The problem is of 
course that these functions are not known. It is 
known, however, how they behave both at low 
and at high values of Q 2 thanks to Chiral Pertur- 
bation Theory and the Operator Product Expan- 
sion, respectively. So, if we can find a reliable way 
to interpolate between the two, our job is done. 
At this point, large-A^ comes to help. Since 
in the large- N c limit both functions Wdg(ll)rr 
are made out of an infinity of zero-width reso- 
nances, they are meromorphic functions. In other 
words, they only have isolated poles, but no cut. 
Furthermore, they are order parameters of chi- 
ral symmetry and would vanish were it not for 
its spontaneous breakdown. It follows, therefore, 
that perturbation theory yields a vanishing result 
to all orders in a s and, in particular, there is no 
parton-model logarithm. With all these consider- 
ations, our choice for the interpolator is then the 
most natural one, namely a meromorphic func- 
tion with a finite number of poles. The position 
of these poles will be identified with the reso- 
nance masses, and the unknown residues will be 
determined so that the interpolator reproduces 
the low- and high-Q 2 expansions given by ChPT 
and the OPE. In mathematics this is known as a 
rational approximant. Therefore, our interpola- 
tor constitutes an approximation to the large- N c 
curve. 

I would like to emphasize that this approxima- 
tion is systematic: the more terms in the OPE 
and ChPT are known, the more resonances can 
be included in the interpolator 4 . Although the 
solution to large- N c QCD is not known, it is plau- 
sible that such an approximation may do a nice 
job. For one thing, we are interpolating a QCD 

4 The interpolator which matches just the first term in 
the OPE is what we have sometimes called the "Minimal 
Hadronic Approximation" 3 because it is the simplest one 
which guarantees the correct short distance properties for 
the matrix element. 



Green's function in the euclidean. Therefore one 
should not expect a lot of "structure"; i.e. reso- 
nances do not show up as "peaks", unlike in the 
minkowski region. For another thing, we interpo- 
late in the gap between the regime governed by 
ChPT and the OPE, so the gap does not seem 
very large! Of course, in the end, one will have 
to judge by the results obtained. 

What are the generic properties which one can 
expect for functions such as Wdg(ll)rr{Q 2 ) ? 
Figure 1 is a schematic view of the expected typi- 
cal profile for one such generic QCD Green's func- 
tion, here called G(Q 2 ). As we can see, at Q 2 = 
the function reaches a value determined by a typ- 
ical chiral parameter, such as Fq or (i/iip), whereas 
at large Q 2 it falls off like an inverse power of Q 2 . 
The turning point is given by a typical resonance 
mass, of the order of 1 GeV. 



Q 2 G(Q 2 ) 




Mr — 1 GeV 2 Q 



Figure 1. Schematic view of the QCD Green's 
function. The interpolation between ChPT and 
the OPE is represented by a (red) dashed curve . 

Let us come back to the matching condition © . 
In trying to build some intuition, it is important 
to keep track of all the different hadronic scales 
involved in the problem since in QCD there is a 
very clear hierarchy, namely 

F < (V^) 1/3 < Mi?. -lGeV, (10) 

and this determines the size of the final contribu- 
tion. For instance, let us now estimate the size 
of the unfactorized contributions relative to the 
factorized ones in the case of gs- Equation JSJ 
shows that the factorized contribution from Q4 is 
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of order unity. In order to estimate its unfactor- 
ized contribution one can use that the shape of 
W LLRR (Q 2 ) is like that of G(Q 2 ) in Fig. 1 and 
therefore 



J dQ 2 Q 2 W LLRR (Q 2 ) ~ F 2 M 2 , 



(11) 



up to logarithmic factors. Consequently the ratio 
of contributions 



unfactorized 
factorized 



16tt 2 F 2 



(12) 



which, although suppressed at large N c like 
0(1/N C ), is a number of order unity. The im- 
portant point is that this effect is very generic, 
and does not depend on the particular form of 
the four-quark operator. For instance, despite 
appearances, the contribution from Qq also obeys 
(|12|) : the factorized contribution is proportional 
to L 5 (^pip) 2 /F§, but the function W DG rr{Q 2 ) 
has the shape of G(Q 2 ) in Fig. 1 multiplying its 
value at the origin by (L 5 - 5L 3 /2)(^) 2 /.F 4 g|. 
This is a consequence of Qe being made up of 
scalar and pseudoscalar densities rather than left 
and right currents. 

One simple reason behind the estimate given 
in l|12l) is that unfactorized contributions know 
about the scale M R ~ 1 GeV whereas factoriza- 
tion is dominated by chiral parameters, like F 
or (ipip), whose typical scales are smaller. Fac- 
torized contributions do not have "access" to the 
scale M R . In this sense it is something similar to 
the opening-up of a "new channel" in a scatter- 
ing process. Therefore, there is no reason to ex- 
pect the next-to-next-to-leading contribution to 
be even bigger, with the consequent breakdown 
of the large- N c expansion. Since the large con- 
tribution is due to the coming into play of a new 
scale, once all the scales have already appear in 
the problem there should be no more surprises. 

There is another (independent) reason why, in 
the particular case of Qe, one could have expected 
large unfactorized corrections: the factorized con- 
tribution depends on Ls(/i). Although the run- 
ning of L5 with \x is a 1/N C effect, in fact, a small 
change in scale such as going from the rho mass 
to 1 GeV already changes the factorized contribu- 
tion by roughly a factor of 2 ! If one takes this as 



a naive estimate for the 1/N C corrections, one also 
concludes that it is not unnatural for these to be 
large. Indeed, it was found that the unfactorized 
contributions were actually larger than the fac- 
torized ones^J. Within the context of a model, 
this was also found in Refs. [H]. 

Using that the Wilson coefficients can be de- 
composed as|2j 

V t * s V td 



Ci{p)= Zi(jj,) +ryi(p), with t= - 

the coupling constant gs in Eq. 
an imaginary part due to r. 
(Imr = -6.0(5) • lfr 4 [S]) 

Imgs — (3 ± 1) X Imr , 



u*v d ' (13) 

"us v ud 

||SJ| picks up 
We obtain||] 

(14) 



where the error has been estimated by varying 
the quark condensate (which is the source of the 
biggest uncertainty), (^) 1 / 3 (2 GeV) = (0.240- 
0.260) GeV . In Ref. it was also obtained that 5 



Im (e g e 



(1.6 ± 0.4) • 10~ 6 GeV 6 x Imr .(15) 



With (fTBTl and (fT4*|) one can get to an estimate for 
e'/e®: 



e[ = lm(V t *V td )G F LJ 
£ 2|e||ReA | 



pW(l-filB)--P (2 



(16) 



where 

p(0,2) . 



E 

i=4,6,8 



Vi it 1 ) { C 71 " 71 ") 0,2 1 Qi (AO I K°) 



(17) 



and fijg = 0.16±0.03 ||| is a term induced by 
isospin breaking. Using the above value for the 
condensate, 1/cj = 22.2 and the physical value for 
RcAo, we obtained in Ref. @] 

E i 

- ~ (2 ±2) x 10~ 3 . (18) 

£ 

However, using the recent determination J7ib = 
0.06 ± 0.08 ^01 an( i including final state interac- 
tions as in Ref. JI] one obtains instead, 

E i 

- ~ (5 ±3) x 10~ 3 . (19) 

£ 

This is to be compared to the experimental num- 
ber, i.e. e'/e = (1.66 ±0.16) x 10~ 3 12 . The les- 
son is that, due to the difference between the two 



3 See Ref. |S| for other analyses. 
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terms in (I16fl , small corrections to either term get 
amplified making this observable extremely diffi- 
cult to control. Matching the level of precision 
of the experimental result is going to be a very 
arduous task. 

At present, the lattice determination of e'/e 
|13| is smaller than the experimental result by 
roughly a factor of three and differs in sign. Our 
previous considerations could be of help in un- 
derstanding this situation. First of all, on the 
lattice one does not really have QCD, but a 
quenched version of it. This quenched version 
has some ghost quarks which, although spin one- 
half particles, are bosons. They are responsible 
for canceling the quark contribution in any inter- 
nal quark loop. The presence of these ghosts dra- 
matically changes the flavor symmetry of the the- 
ory from the usual SU(3)l x SU(3)r to a graded 
SU(3\3) L x ,9r/(3|3) P [I3|. As a consequence, an 
operator like Qq in Eq. J7J has no longer well- 
defined transformation properties under the fla- 
vor group but it can be decomposed as|15| 



Q 6 = 



a QNS 



(20) 



Q QNS 



2(^7X)(9«7^,-9fl7^), 



where q is summed over u, d, s, and q = u,d,s 
are the bosonic (ghost) quarks. Whereas Q§ S 
is a singlet under the quenched group SU(3\3)r, 
Q® NS is not. This fact changes their correspond- 
ing chiral properties completely. To leading order 
in ChPT, these operators may be represented by 



Q QNS 



agistor (AC^Cn 



N 



+ag' 1} str (2B A(UM + MU ] )) , 
/ 2 < s str (ACWl/t) , 

~diag(l, 1,1, -1,-1,-1), (21) 



with qAq = sd, M the quark-mass matrix, and 
str is the so-called supertrace. N represents 
the non-singlet structure of Qg NS . Comparing 
Eqs. and 121fl . one sees that the couplings 
a qi2 nav e a corresponding counterpart in the un- 



( 8 1) 

quenched theory (the weak mass term of a^ 2 ' 
is not written in Eq. (J2J), but a^ s is a total 
quenching artifact. Presently, is neglected 

in lattice simulations. 

Notice that Q§ NS is very similar to the elec- 
troweak penguin Qg, if one makes the replacement 
Q — » N. It follows that must obey the same 



matching condition g ew does, and this leads to[7] 6 

(22) 



.NS 



1 + 



1 



where B = -{ip^/Fg. 

As to the operator Qq S , its contribution to 
ag' 1 ^ is also given by Eq. JSJ in the leading- 
log approximation. One can neglect the contri- 
bution from Q4 in this expression because it is 
numerically small and, furthermore, Qg and Q4 
no longer mix in the quenched theory at leading 
log. One then concludes that the large unfac- 
torized contributions which were dominant in the 
unquenched theory are now vanishing due to the 
fact that rip is effectively zero, as a consequence of 
the quark-ghost cancelation One predicts, 

therefore, that the result is given by the factor- 
ized contribution, which is a much smaller value 
for a ^ than in the unquenched theory : 



a; n ' ' = -I6L5F0 i?o 



(23) 



Moreover, if one compares a^ s 



,,(8,1) 



finds 



with Oiq{ ' OIK' 



a 



I6L5 



60, 



(24) 



where we have used that L5 ~ 10~ 3 . This is not a 
small number and obviously questions the validity 
of neglecting a% s in lattice simulations of e'/e as 
is currently done. 

I would like to conclude by coming back to 
unquenched QCD and turn to the AI = 1/2 
rule. This rule can be expressed by the fact that 
Re g^ xp ~ 3.3 is much larger than g^ p — 0.23, af- 
ter subtraction of the chiral corrections ^7j. Any 
systematic framework must also face the calcula- 
tion of these values. 



3 Ref. 1161 presents another derivation of this result. 
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Although -due to mixing- below the charm 
mass all operators contribute to Re g$, at fi = m c 
only given by 



Q 2 
Qi 



4(sL7 M «i)(tti7/i^i) 
4(sL7 Al di)(ML7 AI iii) 



SGR00188 and by TMR EC-Contracts HPRN- 
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